Phase transitions in the Potts spin glass model 



Giancarlo Franzese 1 ' 2 and Antonio Coniglio 1 ' 2 ' 3 
1 Dipartimento di Scienze Fisiche, Universitd di Napoli, Mostra d'Oltremare Pad. 19 1-80125 Napoli Italy 
2 INFM - unita di Napoli 3 INFN ■ 



00 1 
ON i 
G\ ■ 

5h ! 
Or 

<: 
m : 



43 

o 



c 

o 
o 



> 

m 
o 

o 

oo 
On 
-i— > 



i 

c 

o 
o 



X 



unita di Napoli 

(February 1, 2008) 



sezione di Napoli 



We have studied the Potts spin glass with 2-state Ising 
spins and s-state Potts variables using a cluster Monte Carlo 
dynamics. The model recovers the ±J Ising spin glass (SG) 
for s = 1 and exhibits for all s a SG transition at Tsg(s) and a 
percolation transition at higher temperature T p (s). We have 
shown that for all values of s ^ 1 at T p (s) there is a thermody- 
namical transition in the universality class of a ferromagnetic 
s-state Potts model. The emciency of the cluster dynamics is 
compared with that of standard spin flip dynamics. 

PACS numbers: 75.10.Nr, 64.60.Ak, 02.70.Lq 



I. INTRODUCTION 

In nature there are many examples of glassy systems, 
i.e. complex systems that exhibit a very slow dynamics 
that prevents them to reach the equilibrium state. In this 
class there is a large variety of systems, like real glasses, 
spin glasses, supercooled liquids, polymers, granular ma- 
terial, colloids, ionic conductors, orientational glasses, 
vortex glasses M. A common feature of these systems 
is the frustration, that is a competition due to geometry 
or energy constrains. 

Experiments on glasses Q show that in a temperature 
driven transition precursor phenomena starts at a tem- 
perature well above the ideal glass transition. This tem- 
perature depends on experimental conditions and is the 
onset of some dynamical anomalies, like non-exponential 
relaxation functions, anomalous diffusion, cooling rate 
dependent density. Analogous phenomena are present in 
spin glasses (SG) ||, that are a generalization of Ising 
model where random distributed ferromagnetic and an- 
tiferromagnetic interactions give rise to frustration. In 
particular, experiments [Q and numerical simulations 
pj have shown non-exponential autocorrelation functions 
below a temperature T* well above the transition tem- 
perature TsGi while above T* only exponential relaxation 
functions are seen. 

To explain this phenomenon, Randeria et al. || have 
suggested that in the SG the non-exponential regime 
starts at the Griffiths temperature T c , i.e. the criti- 
cal temperature of the ferromagnetic model, due to the 
presence of randomly large unfrustrated regions, allowed 
by the disorder in the SG interactions. Campbell et 
al. Q have proposed an alternative hypothesis in which 
the onset T* of non-exponential behavior coincides with 
the percolation temperature T p of Fortuin-Kasteleyn ||] 



Coniglio-Klein @ (FK-CK) clusters (defined below) . The 
idea of Ref . j?J is that in the SG the accessible phase space 
is simply connected above T p , while it is not below T p . 
Therefore for T < T p a local Monte Carlo (MC) dynam- 
ics performs a random walk on a ramified percolating- 
like structure with many time scales, giving rise to non- 
exponential relaxation functions. 

Since T p is less than, but close to, T c and T* is difficult 
to localize, it is not possible to exclude the T* = T c 
hypothesis. However numerical results [[W 11 on fully 
frustrated models, without disorder, where the Randeria 
et al. argument does not apply, are consistent with the 
T* = T p scenario. In Rcfs. |uj[ll] this result is given 
both for spin-flip and bond-flip dynamics. The latter is 
strictly related to the bond frustrated percolation (FP) 
where frustrated loops, i.e. closed paths of connected 
bonds covering an odd number of negative interactions, 
have zero weight [JuJ. The FP model and the ±J SG 
model can be recovered as particular cases of a 2s-state 
Potts spin glass (PSG) |12] respectively for s = 1/2 and 
s = 1. 

In this paper we will show numerical results on the 
static properties of the PSG model in 2D. For any s the 
model exhibits a SG transition at a temperature Tsg(s) 
(in 2D Tsg = for any s) and a FK-CK percolation 
transition at higher temperature T p (s). For s =/= 1 the 
higher transition corresponds to a real thermodynamical 
transition of a s-state Potts model Jl3|] . To make the 
algorithm faster we use the Swendsen and Wang (SW) 
(TJ| cluster MC dynamics that prevent the slowing down 
for temperatures near T p (s). Dynamical properties of 
the model for s = 2 are given in Ref. Jig] , where we have 
shown that autocorrelation times diverge at T p and, as 
in SG model, non-exponential relaxation functions are 
present below T p . 

These results on SG model @, FP model (To), fully 
frustrated systems [0 and PSG model suggest that 
the FK-CK percolation may play a role in the context of 
precursor phenomena, since below T p the frustration is 
present at all length scales by means of FK-CK percolat- 
ing cluster, that cannot include frustrated loops. 

In Sec. II we define the Hamiltonian and review some 
theoretical results. In Sec. Ill we present the numerical 
results for the model in finite dimensions. In Sec. IV we 
define the SW dynamics and compere it with the spin flip 
dynamics, verifying the efficiency of the cluster dynamics. 
In Sec. V we give the conclusions. 
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II. HAMILTONIAN FORMALISM 

The Potts spin glass (PSG) model is defined by the 
Hamiltonian 



H 



(i,3) 



-ijSiSj 



I) 



2] 



(1) 



where to each lattice site is associated an Ising spin Si — 
±1 and a s-state Potts spin cr, = l,...,s. The sum 
is extended over all NN sites, e;.j = ±1 is a random 
quenched variable and J is the strength of interaction. 
The model is a superposition of a ferromagnetic s-state 
Potts model and a ± J Ising SG J3| and for 5 aiaj = 1 
(i.e. s = 1) recovers the ±J Ising SG Hamiltonian. 



Following Ref. 
ters on this mode 



it is possible to define FK-CK clus- 
., activating a bond, between NN sites 
with both SG interaction and Potts interaction satisfied, 
with probability 



p = l 



-2sJ/k B T 



(2) 



and defining a cluster as the maximal set of connected 
bonds. 

For a given set of interaction {ejj} it is possible to 
shown (12) that Z can be expressed in terms of bond 
configurations C 



Z{e itj }= £ 



,-H/k B T 



(3) 



{S„a,} 



where W S (C) = if C includes any frustrated loop, oth- 
erwise 



W s (C)=p\ c Hl~p) lAl (2s) N ^ 



(4) 



where p is given in eq. (||), N(C) is the number of clus- 
ters in the configuration C, \C\ is the number of bonds 
and |C| -I- I A I is total number of interactions. Let us ob- 
serve that, while the Hamiltonian (jl]) is defined only for 
integer values of s, the eq. (||) is meaningful for every 
values of s and that for s = 1/2 the eq. (H) gives the par- 
tition function of bond FP where a bond configuration 
without frustrated loops has a weight W(C) = e^ M ' c 
with /3/j. = h\(e^ J — 1), while a bond configuration with 
frustrated loops has a zero weight. Furthermore in the 
limit s — * eq. (|3j) gives the partition function of the 
tree percolation [ |13| ] where any bond configuration with 
a loop is excluded. 

Exact renormalization group (RG) analysis on hierar- 
chical lattice 16 1 has predicted for the PSG two criti- 
cal temperatures. The lower temperature Tsg(s) corre- 
sponds to a SG transition in the universality class of ± J 
Ising SG, and the higher T p (s) to a percolation transition 
in the universality class of a ferromagnetic s-state Potts 
model. Same results are given for the fully frustrated 
version of the model studied with a mean field approach 



Looking at the partition function (Q) one should ex- 
pect a singularity at T p (s), for any s / 0, due to the 
singularity in the number of cluster N(C). Nevertheless 
this singularity has never been observed in the case of SG 
(s = 1). In fact, the RG calculations in the case s = 1 
show a singularity at Tsg for the SG free energy and no 
sing ularity at T p (l). This result is interpreted in Ref. 
|16| supposing that the free energy of the Hamiltonian 
(Jfi) has the form 



F s (T)-F s (T p )~A( S )(T-T p (s)) 



2-a(s) 



where A(s) is an amplitude which vanishes for s 
a(s) is the specific heat exponent. 



(5) 
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III. MONTE CARLO RESULTS 

We have done our simulation using the Swendsen and 
Wang (SW) Monte Carlo cluster dynamics Q described 
in Sec. IV. As we will show, the SW dynamics is faster 
then standard local dynamics, like spin flip dynamics, 
but suffers of slowing down near the SG critical temper- 
ature. Nevertheless, since the SG transition in 2D occurs 
at Tsg = jl8) and we are interested in study the sys- 
tem near the percolation transition at T p > Tsg, the SW 
dynamics is particularly indicated. 

We have performed numerical simulation of the PSG 
model for s = 2, 7, 50 on a 2D square lattice with linear 
sizes ranging from L — 10 to 60 lattice steps and with 
quenched random interaction configurations {fir/}- Defin- 
ing as MC step an update of all the spins of the system, 
we have discarded the data of the first 7 500 MC steps 
and have collected data over 15 000, 25 000 or 50 000 
MC steps, depending on the temperatures and sizes. 

For each s we have calculated the Binder parameter 
for the energy density E defined as 



V = 1 



(E 4 ) 



(G) 



where the symbol (•) stands for the thermal average. This 
quantity allows to localize the transition and to distin- 
guish between first order and second order phase transi- 
tion. In fact, for a second order phase transition in the 
limit L— >ooitisl/ = 2/3 for all temperatures, while 
for a first order phase transition it is 



V 



1 (E+ - E. 
3 _ 



\E + + E_Y 



{E\ + E 



2 \2 



(7) 



where V mm is the minimum value of V (occurring at the 
phase transition temperature) and E + — E_ is the energy 
jump, related to the latent heat, at the same tempera- 
ture. 

To estimate the thermodynamical critical exponents 
for the second order phase transition we have measured 
the Potts order parameter 



2 



M = 



s maxj(M,) — 1 



s- 1 



(8) 



(where i = 1, . . . s, Mi is the density of Potts spins in the 
ith state), the susceptibility 



X 



(M 2 ) - {Mf 



N 



1/2 



(9) 



(where N is the total number of Potts spins) and the 
specific heat 



C H = 



(E 2 ) {Ef 



-,1/2 



N 



(10) 



Furthermore to estimate the percolation critical expo- 
nents we have calculated the percolation probability per 
spin 



P = 1 - kllk 



(11) 



where k is the cluster size and ilk is the density of clusters 
of size k, the mean cluster size 



and the number of clusters 



N c = n k . 

k 



(12) 



(13) 



From standard scaling analysis applied to finite sys- 
tems [BOf we expect for M 



M~L-P!»f M {{T-T s )L 1 l 1 ') 



(18) 



where Jm(x) is an universal function of the dimensionless 
variable x. Analogous scaling functions are expected for 
the other thermodynamical quantities. Tuning the values 
of critical exponents and T s it is possible to verify the 
scaling hypothesis, as eq.([l8|), from the MC data. The 
values for which the data collapse give the estimates of 
the critical exponents and of T s . 

In Figs. H |we show the data collapses for system 
sizes L — 10 -i- 60. The estimated scaling parameters are 
|H given in Tab. | 

The estimated critical exponents are compatible, 
within the errors, with the expected values for a Potts 
model with s = 2 state, i.e. an Ising model, in 2D: a = 0, 
(3=1/8 = 0.125, 7 = 7/4 = 1.75 and v = 1 @. The 
estimated critical temperature is fc_eT s / J = 2.95 ± 0.15. 

For the FK-CK percolation quantities we have de- 
scribed the critical behavior of P and S introducing a 
percolation set of critical exponents (a p , j3 p , 7 P and v p ) 
defined by the relations 



\T-Tp\~ v ' 



(19) 



where £ p is the connectedness length of the clusters and 
T. p = limi^ 00 T p (L) with T p (L) finite size percolation 
temperature, 



P~ \T-Tpf* 



@vl v v 



(20) 



A. Results for s — 2 

In Fig.jl] we show the Binder parameter for the case 
s = 2 for system sizes L = 10 -j- 60. It is possible to 
see that V for small sizes has a minimum at ksT/J ~ 
3.0 and that for grater sizes it becomes constant for all 
temperatures, reveling a second order phase transition. 
Therefore we can make standard scaling analysis p0| for 
the thermodynamical quantities. 

In particular, by definition of critical exponents v it is 

t~\T-T,\- v (14) 

where £ is the correlation length and T s = lhriL^oo T S {L) 
with T S (L) finite size critical temperature of the PSG 
model. Analogously from the definitions of the other 
critical exponents f3, 7 and a we get 

M~ \T-T„f ~t-M v , (15) 
X ~|T-T s r (16) 
C H - \T - T s \- a ~ Z a/U ■ (17) 



S~ \T-T P \-^ -Cp p/,/p • (21) 

Standard scaling analysis for finite systems |23j is applied 
also in this case and the results are summarized in Figs, 
sl ^. The estimated scaling parameters are given in Tab. 

All the estimated exponents are compatible, within the 
errors, to the corresponding thermodynamical parame- 
ters for the 2D Ising model and the numerical estimate for 
the percolation temperature is ksTp/J = 2.925 ± 0.075 
consistent with the estimates of T s . 



B. Results for s — 7 and s = 50 

In Fig. [?] we show the Binder parameter V for the PSG 
model with s = 7 and s = 50 for system sizes L = 10—50. 
The fact that V has a non vanishing minimum for every 
size revels that there is a first order phase transition. In 
this case there is no diverging length, therefore the scaling 
analysis cannot be applied. This kind of transition is 
characterized by the finite size relations pTt 

C H (T S (L),L) ~ max T (C7tf(T,L)) ~ L D (22) 



3 



(where D is the Euclidean di men sion) for the maximum 
of finite size Cr{L) (see Tab. Ill) and the relation 



2max(A) — Tmax(oo) ~ L 



-D 



(23) 



where T m ax(^) is the temperature of the maximum of 
Ch (or of the mean cluster size S) for the size L and 
Tmax(oo) is the corresponding value in the thermody- 
namical limit, i.e. the corresponding transition temper- 
ature T s (s) (or T p (s)). Therefore T s (s) and T p (s) can 
be evaluated by linear fits with one free parameter. The 
data are given in Tab. IV and the results are for s = 7 



T p = 7.5 ±0.1 and for s = 50 T s = T p = 35.0 ±0.1. 



The results are summarized in the phase diagram in 
Fig.|§|: For every s the high temperature phase is disor- 
dered and non percolating; decreasing the temperature 
there is a second or first order phase transition (depend- 
ing on s) at T p corresponding to the percolation of FK- 
CK clusters and to the ordering of Potts variables; at 
lower temperature there is the SG transition (that in 2D 
occurs at T — [jlql). For a fixed realization of {e^} it 
is possible to show pjj that each critical point is charac- 
terized by a diverging critical length. At T p (s) diverges 
the linear size of FK-CK clusters associated with the pair 
connectedness function, while at Tsg(s) diverges the lin- 
ear size of correlated regions. 

It is interesting to note that the behavior of T p (s) can 
be obtained from the exact expression of the transition 
temperature of a ferromagnetic 2s-state Potts model |l3| ] 
only by renormalizing the number of states, i.e. 



(24) 




with a = 0.803 ± 0.003 (choosing J = k B )- 



IV. COMPARISON WITH SPIN FLIP DYNAMICS 

The MC dynamics used to study the equilibrium prop- 
erties of the PSG model is the Swendsen and Wang (SW) 
cluster MC dynamics. The SW dynamics is performed 
in two steps. The first steps is to construct the FK-CK 
cluster configuration C, given an Ising and a Potts spin 
configuration {Si,<Ji\, activating bonds with the proba- 
bility in eq. (pj) between NN sites when both Ising and 
Potts spins satisfy the interaction and with zero proba- 
bility otherwise. The second step consists in reversing all 
the spins in a cluster at the same time with probability 
1/2, for each cluster. The sequence of the first and of the 
second steps applied to the whole system constitutes a 
MC step, that is the chosen unit of time. 

This dynamic completely overcomes the problem of 
critical slowing down for the unfrustrated spin models 
jl4j], while turns out to suffer of diverging correlation 
time if applied to frustrated systems near a critical point. 
This inefficiency is a consequence of the fact that the FK- 
CK clusters used in the SW dynamics do not represent 



anymore in frustrated models the regions of correlated 
spins near a critical point [jl2| and their percolation tem- 
perature T p is greater than the critical temperature. In 
particular, this is true for the SG model for which an ef- 
ficient cluster MC dynamics does not yet exist except for 
2D ]2q] , while efficient cluster dynamics have been pro- 
posed for systems with frustration but without disorder 



Nevertheless in SG for temperatures well above the 
critical temperature T$g and near T p the SW dynamics 
is still efficient, consistently with the general observation 
that the cluster dynamics are efficient at least for tem- 
peratures above the percolation temperature p8| . 

On the other hand in Ref. [jl5[ we have shown that the 
local spin-flip (SF) MC dynamics fo[ in the case of PSG 
model with s = 2 is characterized near T p by diverging 
correlation times. To compare the efficiency of SW dy- 
namics to that of SF dynamics we have studied for the 
case s = 2 the correlation functions at the equilibrium, 
that for a generic observables A is defined as 



fA(t) 



(S(t + t )S(t Q )) 



(25) 



where 8(t) = A(t) — (A) and to is the equilibration time. 
As observables we have choose the Potts order parameter 
M and the energy E of the whole system. 

We have also studied the time dependent nonlinear 
susceptibility for a quenched interaction configuration 



Xsc(t) = — 



(26) 



(where N is the total number of spins). The normalized 
correlation function is 



fx 



Xsajt) - xsajt = oo) 
Xsg(0) - xsc(t = oo) 



(27) 



with X sg(0) = N. 

For both the SW and the SF dynamics we have mea- 
sured the integral correlation time defined as 



Tint 



lim 



1 



E 

t=0 



/(*) 



(28) 



where / is the generic correlation function. We have 
considered systems with lattice sizes L < 30 at tempera- 
tures above and below T p . The data for the SF dynam- 
ics are averaged over 32 different quenched interaction 
configurations, since the local updating of this dynamics 
strongly depends on the local fluctuation of the frustra- 
tion. On the other hand the results on global SW dy- 
namics turns out to be "robust" respect the interaction 
configuration average, in the sense that the fluctuations 
of Ti n t are within the errors estimated on the basis of a 
single interaction configuration analysis. 
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The simulations have been done with an annealing 
method, i.e. with a slow cooling of the system at each 
temperature. For the SW cluster dynamics 5 10 3 MC 
steps turn out to be enough to equilibrate the system at 
the considered temperatures and the averages are done 
using the the data for the following 5 10 4 MC steps. For 
the SF dynamics we have discarded the first 10 4 MC steps 
(defined as the local update of any spin in the system) 
and recorded the data for 5 10 5 MC steps. 

In Tab. [v| we show the results for L — 30. Analogous 
results have been found for smaller systems. 

The data show that, while the SF correlation times 
for M and E grow abruptly near T p where both ther- 
modynamical and percolation transition occur, the SW 
correlation times only show a slow trend to increase for 
decreasing temperatures, being smaller than the corre- 
spondent SF data at least of an order of magnitude. Even 
for t x , that for SF starts to be non-zero below T p , the 
SW dynamics shows smaller correlation times. For tem- 
peratures well below T p it is possible to see that the SW 
dynamics is characterized by long autocorrelation times, 
as the SF dynamics. Therefore, at least for temperatures 
not too much below T p , the SW dynamics turns out to be 
more efficient than local SF dynamics. In particular, near 
the thermodynamical transition at T p the SW dynamics 
completely overcomes the critical slowing down problem 
in the PSG model, as for the unfrustrated models, even if 
frustration is present via the random interactions of the 
Ising spins. 



V. CONCLUSIONS 

Precursor phenomena characterize the paramagnetic 
phase of spin glass. In particular, experiments and local 
Monte Carlo simulations show the presence of stretched 
exponential autocorrelation functions well above the SG 
transition temperature T$g @ ■ The relation of the onset 
T* of these precursor phenomena to any thermodynami- 
cal transition and the localization of T* are still matters 
of debate [f7||l5|l. Many numerical evidences on disor- 
dered Q and deterministic frustrated models [ |lO|JTl| ] have 
shown that T* is consistent with the percolation temper- 
ature T p of the the Kasteleyn - Fortuin and Coniglio - 
Klein clusters. In particular, a generalization of the ± J 
Ising SG to a 2s-state Potts SG, that recovers the SG for 
s = 1, has shown that for s — 2, as for s — 1 the 
T* = T p hypothesis is numerically verified EH]. In this 
paper, using very efficient cluster Monte Carlo dynamics, 
we have shown that in 2D for any s ^ 1 the percolation 
transition corresponds to a real thermodynamical transi- 
tion in the universality class of the s-state ferromagnetic 
Potts model. In particular, we have considered the case 
s = 2 and the cases s = 7, 50 where at T p (s) a second 
order and, respectively, a first order phase transition oc- 
curs. Exact renormalization group calculations on hier- 
archical lattices fill and mean field analysis for a version 



of the model without disorder p7| have shown the same 
scenario. 

All these results suggest that the percolation transition 
may play a role in the precursor phenomena even in the 
SG case (s = 1), where no thermodynamical transition 
occurs at T p . This idea arises from the observation that, 
even in SG, below T p the frustration starts to be man- 
ifested on all the length scales by means of the FK-CK 
clusters, that cannot include frustrated loops. Therefore 
the scenario presented is that for s ^ 1 at the percolation 
temperature T p there is a thermodynamical transition 
with associated dynamical anomalies that " vanishes" for 
s = 1 leaving the dynamical behavior unchanged. 
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TABLE V. Integral correlation times for a PSG model 
with s = 2 and L = 30 for Swendsen and Wang (SW) clus- 
ter MC dynamics and local spin-flip (SF) MC dynamics for 
temperatures above and below T P (L = 30) ~ 2.95J/fc s . 



k B T/J 


2.75 


3.00 


3.25 


tm(SW) 


3.08 ± 0.02 


3.65 ± 0.03 


1.81 ±0.01 


tm(SF) 


76.8 ±0.3 


575.8 ± 0.6 


180.3 ± 0.2 


te(SW) 


9.32 ± 0.02 


8.53 ±0.07 


4.45 ± 0.04 


te(SF) 


19.0 ±0.1 


67.5 ±0.1 


20.9 ±0.5 


T X (SW) 


9.93 ±0.07 


2.423 ± 0.002 


1.585 ±0.001 


r x (SF) 


16.89 ±0.02 


7.06 ± 0.03 


3.97 ±0.01 



TABLE I. Estimated critical exponents and critical tem- 
peratures T a for thermodynamical quantities M, \ an d Ch 
for the PSG model with s = 2. 





v a/v 




7/1/ 


k B T s /J 


M 


0.9 ±0.3 


0.13 ±0.06 




3.0 ±0.1 


X 


1.0 ±0.3 




1.75 ±0.10 


2.9 ±0.1 


C H 


1.0 ±0.5 0.0 ±0.1 






2.9 ±0.1 



TABLE II. Estimated percolation exponents and percola- 
tion temperature T p for P and S for the PSG model with 
s = 2. 





Up 


Pp/vp 




k B T p /J 


P 


0.9 ±0.2 


0.10 ±0.03 




2.95 ±0.05 


s 


0.95 ±0.15 




1.6 ±0.2 


2.90 ±0.05 



TABLE III. Maxima of Ch for s = 7 and s = 50 for 
L = 30, 40, 50. 



L 


30 


40 


50 


max-rCtf (s 


= 7)/L a 


27 ±3 


26 ±3 


33 ±3 


maxrCii(s 


= 50)/L 2 


7.4 ±0.6 


8.4 ±0.9 


7.3 ±0.7 



TABLE IV. Temperatures (in J/k B units) of maxima of 
Ch and S for s = 7 and s = 50 for L = 10, 20, 30, 40, 50. 



L 


10 


20 


30 


40 


50 


Tmax c H ( s 


= 7) 


7.8 


7.6 


7.5 


7.5 


7.5 


Tmax s(s — 


7) 


7.9 


7.7 


7.5 


7.5 


7.52 


Tmax c H ( s 


= 50) 


36.56 


35.57 


34.92 


34.90 


35.00 


2max s(s = 


50) 


36.56 


35.60 


35.00 


35.00 


35.00 
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FIG. 1. PSG model for s — 2: Binder parameter V vs. 
dimensionless temperature fc^T/J for L = 10 -j- 60. Errors 
are smaller than the symbol sizes. Lines are only guides for 
the eyes. 
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FIG. 3. PSG model with s = 2: Data collapse for x f° r 
systems sizes L = 10 -j- 60. Temperatures ate in J/fcs units. 
The scaling parameters are given in the figure. 
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FIG. 2. PSG model with s = 2: Data collapse for M for 
systems sizes L = 10 -j- 60. Temperatures ate in J/fcs units. 
The scaling parameters are given in the figure. 
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FIG. 4. PSG model with s = 2: Data collapse for C* H for 
systems sizes L = 10 -t- 60. Temperatures ate in J/fcs units. 
The scaling parameters are given in the figure. 



7 



-/?/i/=-0.1 ±0.03 

i/=0.9±0.2 

Tp = 2.95±0.05 




T 

o 
□ 

A 


o 



L-- 
L-- 
L-- 
L-- 
L= 
L= 



:60 

^30 

:20 

= 10 



200 400 ,600 
(T-Tp)L (1/ ^ 

FIG. 5. PSG model with s = 2: Data collapses for P 
for systems sizes L = 10 -j- 60. Temperatures are in J/fcs 
units. The scaling parameters are given in the figure. In 
figure indexes p are omitted for the critical exponents. 
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FIG. 7. PSG model: Binder parameter V vs. dimen- 
sionless temperature fcsT/J for (a) s = 7 and (b) s — 50, 
for L = 10-7-50. Insets show the particulars for L — 40, 
50. Where not shown, the errors are smaller than the symbol 
sizes. Lines are only guides for the eyes. 
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FIG. 6. PSG model with s — 2: Data collapses for S 
for systems sizes L = 10 -j- 60. Temperatures are in J/fc_g 
units. The scaling parameters are given in the figure. In 
figure indexes p are omitted for the critical exponents. 
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FIG. 8. PSG model: Numerical phase diagram in 2D. 
The data are fitted with T p /(2s) = a/(ln(l + \/2sa)) with 
a = 0.803 ± 0.003 (choosing J = k B ). Data for s = 1/2 and 
s = 1 are from Ref. [25]. Where not shown, the errors are 
smaller than the symbol sizes. 
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